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INTERACTION  OF  DIFFUSION  AND  BOUNDARY  CONDITIONS 
by  Jack  K..  Hale  and  Carlos  Rocha 

ABSTRACT 

For  systems  of  reaction  diffusions,  the  existence  and  behavior  of  the 
solutions  on  the  compact  attractor  is  discussed  for  large  diffusion  coefficients  and 
boundary  conditions  which  can  vary  from  Neumann  to  Dirichlet  conditions 


1.  Introduction 


Consider  the  system  of  parabolic  partial  differential  equations  (PDE) 

(1.1)  9u/9t  ■  DAu  +  f(u),  xGfl 

(1.2)  D0u/0n  +  0E(x)u  - 0,  xG0O 
* 

where  u  G  RN,  fl  c  Rn,  n  <  3,  is  a  bounded  open  set  with  dCi  smooth,  D  = 
diag(dlv..,dN),  E  «  diag(ej,...,eN),  each  dj  >  0  is  constant,  e^  00  -*  R  is  continuous,  e ^ 
>  0,  j  =  l,2,...,n,  and  0G(O,*»)  is  constant.  The  function  f:  RN-*RN  is  supposed  to 
be  a  C1,1-f  unction;  that  is,  continuous  and  has  a  Lipschitz  continuous  first 
derivative. 

An  interesting  problem  is  the  following  one:  for  fixed  functions  (f,E),  discuss 
how  the  flow  defined  by  (1.1),  (1.2)  depends  upon  the  parameters  (D,0).  In  a  vague 
sense,  the  (D,0)-space  should  consist  of  two  distinct  types  of  points  -  those  for  which 
the  basic  structure  of  the  flow  does  not  change  significantly  when  one  makes  a  small 
change  in  (D,0)  (the  structurally  stable  points)  and  those  points  for  which  a  small 
change  leads  to  a  change  in  the  basic  structure  of  the  flow  (the  bifurcation  points). 

The  purpose  of  this  paper  is  to  make  a  modest  contribution  to  understanding  some 
parts  of  this  problem.  More  specifically,  we  shall  give  some  conditions  on  (f,E) 
which  will  ensure  that  there  is  a  d0  >  0  such  that,  for  any  d  >  d0,  d  ■  min(dj,  j  - 
1,2,...,N)  and  any  0  [0,*»),  there  is  a  compact  attractor  Bd  q  of  (1.1),  (1.2)  which  is 
upper  semicontinuous  in  D,0  uniformly  for  d  >  d^  0  >  0.  Furthermore,  Bd  q  is  a 
singleton  for  0  >  0q,  0o  sufficiently  large  and  converges  to  an  attractor  for  the 
Dirichlet  problem  for  (1.1).  These  results  complement  the  ones  obtained  previously  in 


the  paper  of  Hale  and  Rocha  [7],  in  which  they  proved  the  existence  and  upper 
semicontinuity  of  BD  0  for  d  >  d0  and  9  in  a  compact  set.  The  new  contribution 
is  the  uniformity  in  0  >  0,  which  permits  one  to  go  from  Neumann  boundary 
conditions  to  Dirichlet  boundary  conditions  for  any  d  >  dQ. 

For  a  scalar  equation  in  one-dimension  and  d  >  d0,  the  types  of  bifurcations  that 
occur  as  one  goes  from  Neumann  to  Dirichlet  conditions  is  also  discussed.  There  is 
some  overlap  in  this  example  with  the  work  of  Conley  and  Smoller  [2  ]. 

The  second  aspect  of  the  paper  deals  with  the  classification  of  points  in 
(D,0)-space  as  structurally  stable  or  bifurcation  points.  In  this  case,  we  attempt  in 
Section  3  a  classification  for  a  scalar  one-dimensional  equation  with  f  a  cubic. 
These  results  overlap  the  ones  of  Gardner  [5]  in  a  special  case.  The  proof  of  the 
classification  relies  heavily  upon  the  transversality  theory  of  Henry  [11]. 

To  describe  the  abstract  results  more  precisely,  we  need  some  terminology.  Let 
X  *  L2(n,RN)  and  define  the  operator  A  *  Ap  g  :  D(A)  -  X  by  A$  -  DA0,  where 

D(A)  -  {u€EWs,j(A,Rn)  :  u  satisfies  the  boundary  conditions  of  (1.2)}. 

Then  A  is  a  sectorial  operator  and  one  can  define  the  fractional  powers  Aa  of  A, 
0  <  a,  and  the  space  Xa  ■  D(Aa)  with  the  graph  norm.  If  n  -  2,  or  3,  n/4  <  a  <  1, 
then  Xa  C  W1,J(n,RN)  n  L*(n,RN)  with  continuous  inclusion.  If  n«  l,a*  1/2,  then 
Xt/2  „  w1,l(fl,R)  n  C(fl,R).  We  assume  below  that  a  is  always  chosen  in  this  way. 
One  can  then  show  that  (1.1),  (1.2)  defines  a  local  C1’1  semigroup  TD  0(t)  on  X“ 
(see,  for  example,  Henry  [10,p.75]). 

For  any  set  B  C  X®,  the  w-limit  set  u(B)  is  defined  as  u(B)  -  r\>0  Cl  Ut>T 
Td  g(t)B.  A  set  B  c  X“  is  said  to  be  invariant  if,  for  any  4>  B,  one  can  define 
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TDg(t)*  for  t£R  and  TDg(t)$E  B  for  tE  R.  A  set  AC  X°  is  said  to  be  a 


compact  attractor  for  (1.1),  (1.2)  if  A  is  compact,  invariant  and  there  is  a 


neighborhood  B  of  A  such  that  u(B)  A. 


Let  Xj(D,0)  be  the  first  eigenvalues  of  -djA  with  boundary  conditions 
dj9u/9n  +  0e-u  =  0,  let  0j(D,0)  be  the  corresponding  unit  eigenfunctions,  j  -  1,2,...,N, 
and  $j(D,0)  be  the  N-dimensional  column  vector  with  $.(D,0)  in  the  jth  place  and 


zero  otherwise  and  let 


AD  e  =  diag(X1(D,0) 

»•••»  ^  (D,0», 


•b.O  “  ($j  (D,0),...,$n(D,0). 


The  following  ordinary  differential  equation  (ODE),  corresponding  to  the 


Galerkin  approximation  obtained  by  projection  of  (1.1),  (1.2)  onto  the  N-dimensional 


subspacc  U  spanned  by  the  elements  of  $D  g,  plays  a  fundamental  role: 


(1.3)  dv/dt  -  -AD  g  v  +  f  $D  0  f(»D  0v) 

’  Jn  ' 


For  any  set  B  in  RN,  we  let  Bu  ■  {^d  gv:  v  E  BN).  For  any  sets  B,  C 


Xa,  we  let 


6(B,C)-  sup  dist(x,C) 


A  function  g(X)  from  X  C  Rk  to  subsets  of  X®  is  said  to  be  upper 
semicontinuous  at  XQ  if  lira^^  ^(X^Xq))  «  0.  For  any  set  B  C  X“,  e  >  0,  let 
N(€,B)  be  the  € -neighborhood  of  B.  Let  X®  designate  the  fractional  power  space 
obtained  by  taking  Dirichlet  conditions  for  (1.1)  and  N0(«,B)  be  the  c -neighborhood 
of  a  set  BC  X®.  Our  prinicipal  result  is 

Theorem  1.1.  Lfii  d  -  min(dlv..,dN).  Suppose  there  is  a  compact  set  KCRN  and 
P.QS,ilivC  constants  d0>  0,  60  >  0  such  that  the  ODE  (1.31  has  a  compact  attractor  AD  e 
CK  and  io(N(60,^d  g))  CAD  0  for  each  dad^OaO.  Then,  f or  any  0  <  6X  <  60  and 
an*  <  >  0.  there  is  a  do>do>0  and  a  compact  set  K 1  c  L”  such  that  (1.1).  (1.2)  has 
a  compact  attractor  eCK*  nN(€,Ag0),w(N(81,Ag0)CBD  0  fpr  d)d„,0>O.  Ik 
anragtor  #De  is  Upper  semicontinuous  in  D,  6.  Also  there  is  a  0O  >  0  such  that 
bd  e  is  a  sinRigtpn  0D  e  fci  d  a  d^e  >  e0, 4>D 0  -  4>Dr  as  0  -  •,  0D  •  is  a  solution 
of  the  PirigbJct  problem  fgr  (l.i)  aM  „)) -  (^ 

This  theorem  is  proved  in  Section  2. 

It  is  worthwhile  to  discuss  the  ideas  that  are  needed  to  verify  the  hypothesis  in 
Theorem  1.1.  Suppose  firstly  that  N  -  1,  n  -  (0,1),  and  f(u)  is  a  polynomial  of  degree 
2p+I  with  uf(u)-*-“as|u|-,,*;that  is, 

f(u)-b00uJP+1  +  b1u*P  +  ...  +  bJp+1 

with  b0  <  0.  For  N  -  1,  the  boundary  condition  (1.2)  is  equivalent  to 
dux  -  60ou  -  0  at  x  ■  0 

dux-  03xu  -  0  at  x  -  1 

with  B0  >  0,  Bx  >  0.  The  first  eigenvalue  X  of  -dd2/dx2  with  these  boundary  conditions 
satisfies  0  <  X  <  dn2  and  the  corresponding  eigenfunction  $  -  $(d,0)  can  be  taken  to  be 


positive.  Thus, 


-  def  f1  jp+1  r1 

f(v)  -  #x)f(#x)v)dx  -  E  b,p+H  (  ^+1(x)dx)uj 

*  0  j=0  *0 

so  that  the  signs  of  the  coefficients  in  this  polynomial  are  the  same  as  the  ones  for  f. 

Also,  there  are  constants  k  >  0,  6  >  0  such  that  fV +1  <  k,  j  -  0,1,...,  2p+l  and 

*0 

f  $2p+2  >  8  for  d  >  0,  6  >  0.  So  (1.3)  becomes 
*  o 

(1.4)  dv/dt «  -0Xv  +  f(v). 

Since  0Xj  >  0  and  f(v)/vs  -*  (^  $2p+2)b0  <  8b0  <  0,  it  follows  that  (1.4)  has  a 
compact  attractor  for  every  d  >  0,  0  >  0  and  these  attractors  lie  in  a 
compact  set.  The  hypotheses  of  Theorem  1.1  are  satisfied.  It  is  clear  that  similar 
conclusions  could  be  drawn  for  a  more  general  f  if  the  behavior  of  f  is  appropriate. 

If  N  >  1,  the  hypotheses  are  not  as  easy  to  verify.  For  0=  (0,1),  equations  (1.3)  in 
component  form  arc  given  as 

dvj/dt  -  -0XjVj  +  f  tyxK^xJv . . $n<x)vn)  dx,  j  -  1,2,...,N 

J  o 

Even  though  all  are  positive  on  (0,1),  one  must  assess  their  relative  contributions  to 
the  behavior  of  the  flow  near  v  ■  •  If  all  diffusion  coefficients  are  equal,  then  Xt  =  ... 
■  *N,  ■  •••  *  and  thc  situation  is  much  simpler.  Although  this  topic  clearly  needs 

to  be  investigated  in  more  detail,  it  will  not  be  pursued  in  this  paper. 


2.  Proof  of  Theorem  1.1. 


For  notational  convenience,  we  take  N  =  1,  pointing  out  in  the  appropriate 
places  the  changes  that  are  needed  for  N  >  1.  Also,  let  us  first  assume  n  =  1,  fl  =  (0,1 ), 
so  that  (1.1),  (1.2)  become 


(2.1) 

ut  =  duxx  +  f(u) 

0  <  x  <  1 

(2.2) 

dux  -  680u  =  0 

at  x  =  0 

dux+  001u  =  0 

at  x=  1 

where 

0O,  are  given 

positive  constants  and  0  E  [0,*).  If  H2  *  W2,2(fl,R), 

H1  =  Wx-2(n,R), 

D(Ad  0)  -  {uEH2:  u  satisfies(2.2)) 

Ad,S  "  'du**’ 

then  Ad  0  can  be  extended  as  a  self  ad  joint  operator  in  H1 

f  (Ad  8u)v  -  d  f  ux  vx  +  6[0ou(O)v(O)  +  B^OMl)] 

’  Jo 

defined  for  every  u,vG  H1.  Now  if  we  consider  the  fractional  power  space  X1/2 
defined  as 

X1/2  «  D(Ad  0  +  I)1/2 


tj 


with  the  graph  norm  (see  Henry  [10,  pg.  29]),  we  have: 


'm 


;«g 


m 


|  |  <Ad,0  +  I)1/2u  |  £  =  |  u(Ad  e  +  1)U  =  |  ( Ad  eu)u  +  J  u2 
-d  [  u2  +  f1  u*  +  0Ifiou*(O)  +  31u*(l)] 

•  n  J  n 


Using  the  Sobolev  inequality  u2(x)  <  k  |  |u|  |h1,  we  have 


K,||"||„,‘||(Ad.e*I)'/Iu||Ll<M||u||H, 


where  kj  -  min(d,l),  M  =  [d  +  0(BO  +  61)k]1  ^2.  Since  X  «  L2  and 


I  |(Ad.e  +  i)1,,“|  L  ■  |  |u|Li/i 


we  conclude  that  X1/,2=H1  independently  of  0  (Henry  [10,  p.  167,  exercise  10]). 


Notice  though  that  the  constant  M  in  the  norm  equivalence  grows  with  0,  being 


unbounded. 


We  will  now  consider  the  eigenvalues  and  eigenfunctions  of  Ad  q  and  establish 
uniform  estimates  for  the  eigenvalues.  If  we  let  Ad  q#  =  X$  and  X  =  dy2,  then  the 


eigenfunctions  and  eigenvalues  are 


\  (x)  -  d7j  cos  7.  x  +  06osin  7}  x 


where  yi  are  the  positive  solutions  of 


cotg  7  -  G(7/z) 


g(s)  -  k(s  ■  s*1),  k  -  z  -  ed-He0 


G (y/2) 


Figure  1 


From  this,  we  immediately  obtain  the  estimates 


Xx  e  (0,dnJ);  Xj  >  dir*  for  j  >  2. 


Also,  from  y^tX/d)1/2  and  lim  d'1^,cot(X/d)1/2  -  X'1/*,  we  have 


lim  Xx  -  e(B0+Bx) 

d"**0 


-9 


Next,  we  estimate  |?2  -  yl  |.  If  z  *  n,  then 


d  G<y/n) 

d7 


2k  <L 
n  n 


since  2k  <  1. 


n  b 


3  TT  /  2 


Figure  2 


This  implies  7*  *  7X  >  b  where  b  *  n3/2(n+l)  >  l  when  z  -  n  (see  Fig.  2).  By 
continuity,  there  are  €  >  0,  6  >  0  such  that  7,  -  yx  >  8  for  z6[ff-t,n+t].  If  |z-n|  >  «,  then 
7j  -  7i  >  €.  Therefore,  there  is  a  constant  c  >  0  such  that  7j  *  7X  >  cI/l  for  all  z.  This 
implies  X2  -  >  dc  for  all  d  >  0,  0£[O,*»). 

Now,  as  in  Hale  [8]  or  Hale  and  Rocha  [7],  we  consider  the  decomposition 
■  Y  ®  Y"K  where  Y  ■  span  and  let  T(t)  denote  the  semigroup  generated 
by  Ad  0|Y-*-.  Equations  (2.1),  (2.2)  can  be  rewritten  as 
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v  =  -XjV  +  f(v)  +  P(v,w) 

(2-3)  t 

w(t, •)  =  T  t)w0  +  T(t-s)  Q(v(s),w(s, •  ))ds 


where 


(2.4) 


P(v,w)  «  f  +  w)  -  fXv^)] 

Jo 

Q(v,w)  -  f  (v^j  +  w)  -  ^  $x  f (v$j  +  w) 
f(v)=  f  ^f(^v). 

*  n 


For  d  >  d0,  6  >  80,  the  assumptions  of  the  theorem  imply  that  the  ODE 
(2.5)  v--X1v  +  f(v) 

has  a  compact  attractor  Ad  g  and  it  attracts  a  80-neighborhood  N(8o,/4d0)of  Aaq. 
This  implies  there  is  a  positively  invariant  open  interval  V  containing  Ad  q  and 
Ad  q  attracts  V.  In  the  case  N  >  1,  one  uses  the  converse  theorems  of  Liapunov  as  in 
[8]  to  obtain  a  positive  invariant  open  set  containing  Ad  g  and  which  is  attracted  by 

Ad,6- 

For  each  fixed  8  and  d  >  d0,  it  was  shown  in  [7]  that  equation  (2.3)  has  a  local 
integral  manifold  w«h(v,d,8)  in  a  neighborhood  of  Ad  g. 

From  the  uniform  estimates  of  the  eigenvalue  XJf wc  have,  for  wEY^-, 
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|T(t)wlxi/i  <  k '  t  '1/2  e'dct  |w|x 

|T(t>wixl/,  <  k'cdctiwix 

where  k'  is  independent  of  d,  0.  The  proof  in  [7]  then  shows  the  existence  of  dj  * 
d0  such  that  BdgC  N(«,/4d  g),  u(N(6ltj4d  g))  CBag  for  d  *  dp  0  >  0  provided  that 
we  know  |flL*»  <  k|^xi/2  for  any  0GX1/2,  where  k  is  independent  of  0.  Next,  we 
establish  the  continuous  inclusion  X^CL*  uniformly  in  0  following  Henry  [10]. 
From  the  Nirenberg-Gagliardo  inequality  (Henry  [10,pg.  37]), 

llu||,o  <  C  ||u||®a  Hull2*® 

for  0>  1/4,  and  by  exercise  11,  page  28  of  Henry  [10],  we  have  XaCC°(a>6) 
continuously  if 

Hull  0  <  Cj||(Ad0  + I)u||®  Mull1-®. 

C  '1/1/ 

Thus,  we  need  the  following  estimate,  uniform  in  0: 

||u||  ,  <  K||(Ade  +  I)u|| 

H  L 

If  g  *  (Ad  g  +  I)u,  then  we  can  compute  explicitly  u  as  a  function  of  g.  In  fact, 

duw‘u  “ 


implies 


a 


; 


& 


'5 

iVj 


Hi 


u(x)  -  u(0)  ch  -4  +  u  (O)Vd  sh  -4.  -  -4.  f  g(s)  sh  44  ds 

✓d  x  ✓d  Vd  J0  Vd 

u  (x)  -  u(0)  -4.  sh  -4-+  u  (0)  ch  -4-  -4  f  g(s)  ch  4?-  ds 

x  Vd  Vd  x  Vd  d  Jn  Vd 


From  the  boundary  conditions. 


dux(0)  -  6  B0  u(0),  dux(l)--0BlU(l), 


we  have 


u(l)  -  u(0)  [  ch  -4.  +  0  B0  -4.  sh  -4.1  -  -4.  f  g(s)  sh  45-  ds 

L  Vd  °vd  Vdi  Vdi0  Vd 

ux(l)-u(0)-L[sh-l-_r+  0Bo-4l  ch  J=r]  -  ^  f  g(s)ch45ds. 
✓d  L  Vd  0  ✓d  ✓d  J  d  Vd 


Hence, 


u(0)  ■  [ch  -L.  0(BO+  Bj)  +  ^-sh  ^r<d  +0  *B0B1)|j"1  £  8(s)[ch^ 

+  0B.-ksh  4r^-]ds 
1  Vd  Vd  J 


u(x)' 


»a 


va 


✓d 


(d+e^oBj)  4  sh  4+  e($0+e1)ch  -L 

yd  yd  yd 


h-4*  +60 

Vd 


ds 


-  4  f  g(s)sh  ^  ds. 
✓d  J0  ✓d 


So  finally  we  obtain  ||u]|  ,  <  R(6)  ||g||  ,  where  R(8)  is  a  rational  function  of  6 

L  L 

such  that  for  some  constant  K,  R(0)  <  K.  for  every  0  >  0.  Also  lldu^ll  2  “l|u-g||  2 

L  L 

<  ||u||  }  +  ||g||  2  and  we  easily  obtain: 

L  L 


Hull  J  <  Kllgllj  -  K.ll(Adfl  +  I)u||  ,  . 

H  L  '  L 

This  gives  the  embedding  of  X1^  into  L"  uniform  in  6. 

Our  next  objective  is  to  show  that  Bd  g  is  a  singleton  if  d  >  d^r),  0  >  0o(r) 
with  d^r),  0o(r)  sufficiently  large.  To  do  this,  one  uses  the  following  fact:  for  any 
r  >  0,  there  are  dj  >  0,  0O  >  0  such  that  X^d.O)  >  r  for  all  d  >  dr  0  >  0Q.  Since 
i4d  g  c  K,  a  compact  set  for  all  d^d^OX)  and  i4dg  attracts  N^^g),  there  are 
constants  klt  k2  such  that 

|f(v)|  <  klt  |f '(v)|  <  k,  for  v  Ud>dQ  e>0  y+(N(6o,/<d  0)) 
where  "1+  designates  the  positive  orbit.  If  vv  v2  be  two  solutions  of  (2.5),  then  z- 


Vj  -  Vj  satisfies 
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for  some  £(t)  and  iHvj+OI  <  k3.  Thus,  if  X,(d,0)  >  r  >  ka  for  d  >  dx(r),  0  > 
0o(r),  then  z(t)  -*  0  as  t  -*  •  exponentially.  This  implies  that  the  solutions  of  (2.5) 
approach  an  equilibrium  vo(d,0)  as  t  By  the  same  type  of  argument,  the 
solution  vo(d,0)  is  hyperbolic.  Thus,  Ad  g  is  a  singleton  {vo(d,0»  and  it  is 
hyperbolic  with  exponent  >  r-k2.  This  implies  that  Bdg  is  a  singleton  {4>d  g)  which 
is  a  hyperbolic  equilibrium  point  and  attracts  N(6x,0dg). 

There  is  a  constant  ks  such  that  |f(u)|  <  ks,  for  u€Ud>d^r)  ^  >g  (r)7+N(8x,0d  g). 
Since  4>d  g  is  an  equilibrium  point,  it  follows  that  d2<l>dg(x)/dx1  is  uniformly 
bounded  by  k3.  Thus,  the  set  {<|>d  g,d  >  dx(r),  0  >  0o(r)}  belongs  to  a  compact  set  Kr. 
Let  0j  -  •  as  j  •  be  a  sequence  so  that  ^a,*0  as  j  -»  •  Then,  0d is  an 

equilibrium  solution  of  (2.1)  satisfying  the  Dirichlet  boundary  conditions.  This 
equilibrium  is  hyperbolic  and  therefore  attracts  a  neighborhood  of  itself 
exponentially.  This  neighborhood  can  be  chosen  in  such  a  way  as  to  attract  every  limit 
point  of  {4>d  g.dadjfr),  0  >  0o(r)).  But  this  will  imply  there  is  only  one  limit  point  <l>d  . 
and  completes  the  proof  of  the  theorem  for  N  »  1,  n  ■  1. 

For  N  >  1,  n  -  1,  the  last  part  of  the  proof  follows  in  essentially  the  same  way 
since  one  can  construct  a  quadratic  Liapunov  function  for  (2.6). 

The  case  n  -  2,3  and  arbitrary  N  follows  along  the  same  lines  as  before.  One 
must  obtain  good  estimates  on  the  first  and  second  eigenvalues  of  -A  and  the 
embedding  X®  C  L*  must  be  uniform  in  0.  Because  of  this  last  fact,  all  solutions  in 
the  attractor  can  be  considered  in  L  .  for  0  >  Bo(r)>  there  will  be  a  compact  set  KrC 
L"  which  contains  the  set  {0d  g,  d  *  dx(r),  0  *  0o(r)}  and  ^d>g  "*  ^d«  as  0-*-.  The 
function  0d.  will  satisfy  the  Dirichlet  problem  and  regularity  theory  implies  it  is  in 
x£(n,RN).  Therefore,  we  only  discuss  the  second  eigenvalue  and  the  uniform 
embedding  of  X“  into  L*.  We  also  only  give  the  proof  for  n  -  3  since  obvious 
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changes  will  give  a  proof  for  n  -  2. 

First,  we  establish  uniform  bounds  in  0  for  the  eigenvalues  of  -dA(+BC) 
where  the  boundary  conditions  (BC)  are  9u/dn  +  0eu  -  0  in  dCl  As  in  Hale  and 
Rocha  [7],  we  consider  the  minimum  characterization  of  the  first  eigenvalue 


(2.7) 


min 


jd  J^|?u|j  +  8  |  eu1 


from  which  we  obtain  that  0  <  <  0  |0|-1  e. 

Jan 


If  Xj  =  P  tfj,  then 


=  ^(d/ej-in)'1  f  e  as  d/0-«. 

Jan 

If  Xx  «  dVj,  then 

v1-vl(8/d)-Vio>0  as  0/d-«. 

where  v10  is  the  first  eigenvalue  of  -A  with  Dirichlet  boundary  conditions.  Thus, 
for  any  r  >  0,  there  are  d0  ■  d0(r)  >  0,  80  ■  0o(r)  >  0  such  that  X1(d,0)>rfor  d>  d0, 

0>0O. 

To  estimate  X2,  let  xa«dn,  and  obtain 

M,  =  min([  |7u|J  +  ^-  f  euJ :  f  u*-l,  [  u$.-o) 

d  j an  Jn  Jn  J 


where  ^  is  the  eigenfunction  corresponding  to  Xr  Then  there  exist  positive 
constants  d0,  u  independent  of  0  such  that  >  M  for  every  d  >  dg.  To  see  this, 
consider  n2  ■  Mj(d,0)  and  assume  the  existence  of  sequences  djt  0j  such  that 


djXljj  and  M2(dji®j)  "*  0-  Since  0  <f  l^l2  <  M2»  we  have  f  iv^i2  -0.  Moreover, 

Jn  Jn  - 

denoting  by  $2  the  eigenfunction  corresponding  to  X2,  we  also  have  that  |V$2|J-0. 

f  f 

Then,  from  *  1,  j  -  1,2,  we  have  that  -•  |0|_1/2  contradicting  «  0. 

-*0 

Hence  for  d>d0  we  have  the  estimate  X2>dp  uniformly  in  0. 


We  now  prove  the  following. 


d0  >  °.  i 
X^e.d) 


L  Suppose  n  =  3,  a  >  3/4,  X  =  LJ  aM  X“  =  X“(0,d)  is  the 
r  space  associated  with  -DA  with  boundary  conditions  (1.2).  For  any 
a  constant  k(d0)  such  that  for  any  d  >  d0,  0  E  [0,«)  and  any  u  E 


Proof:  As  for  the  case  n  *  1,  the  essential  step  is  to  consider  the  following  problem: 


Au  =  g  for  x  E  ncR* 


ft 

I 

g 


Bqu  *  0  for  x  E  dd, 


u«  i 

where  B0u  -=  ou/on  +  0u,  0  E  [(),•)  and  to  prove  the  following  uniform 

regularity  estimate: 


(2.8)  ||u||  2  <  M(||g||  2  +  ||u||  ,) 

H  L  L 

where  M  is  a  constant  independent  of  0,  and  then  use  the  fact  that,  for  0  >  0, 
Hull  2  4  c  ||g||  2 ,  (Friedman  [4,  pg.  76]). 

1/  L 

Since  regularity  is  a  local  property,  we  let  E »  1  be  a  partition  of  unity 


subordinate  to  a  neighborhood  covering  of  fi.  Then 


(2.9)  ||u!|*  -  ||  E  p^l*  <  k  E  llpjull*  . 

H  H  H 

If  pjU  has  support  in  the  interior  of  Q,  then  (M  Schechter  [13],  Lemma  7) 

llpjull’  <  CdlPjUll*  +HAp,u|*  A 
H  L  L 

Since  Ap;  u  ■  p4  Au  +  derivatives  of  u  of  order  <  1,  we  have 
||PPiu||*  <2||piAu||J,  +  2c1HuH* 

L  L  H 

and  from  the  inequality  (L.  Nirenberg  [12],  appendix): 

(2.10)  ||u||*  <  €'||u||*J  +  k1(€')||u||J  ,, 

H  H  L 

valid  for  every  €  '  0,  we  finally  obtain 

llPjUll*  ,  <  CdlpjUll*  +  2||PjAuH*j  +  2c  '||u||*  ,  +  2  c  j  1llu||* 

H  L  L  H  L 

<  2C  [||AuH*a  +  (1/2  +  CjkjJHull*^  +  c  'c^ull*^] 

(2.1 1)  HPjuH*  <  K(||g||*  +  ||u||*  +  €||u||*  ). 

n  L  L  H 
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Now,  if  the  support  of  pi  u  contains  a  piece  of  3fl,  we  consider  a  transformation  of 
variables  straightening  up  the  boundary.  We  let  dCL  denote  the  piece  of  dn 

def 

contained  in  the  support  of  v  *  pj u  and  assume  without  loss  of  generality  that  3ft  is 
connected.  Let  0  :  Rn  Rn  be  a  smooth  local  change  of  coordinates  mapping  the 
support  of  v  into  a  ball  Bj  centered  at  the  origin,  and  3ft.  into  B;  n  T,  where 
T  denotes  the  hyperplane  T  ■  {y  ■  (ylt...,yn)  £  Rn  |yn  -  0).  Under  this  local  change 
of  coordinates,  the  initial  problem  is  transformed  into  the  following: 

Lu  =  g  for  y  £  Bj  O  {yn  >  0} 

Bgu*0  for  y  £  BjO  T, 

where  Bgu  ««  (1-0)  3u/3yn  -  0u,  and  L  is  a  linear  second  order  strongly  elliptic 
operator  with  variable  coefficients.  Let  us  denote  by  L„  the  homogeneous 
operator  with  constant  coefficients  which  equals  the  principal  part  of  L  at  the 
origin.  Then,  as  in  M.  Schechtcr  [14],  (proof  of  Lemma  12),  we  may  assume  that  the 
change  of  coordinates  0  (after  a  rotation)  has  the  form  y^  ■  Xj,  j  -  l,...,n-l,  yn  = 
0(xlv..,xn)  such  that,  at  the  point  x0  «  ^(O),  we  have  0(xo)  ■  0  and  also  30/3xn  = 
1,  hence  preserving  Lebesgue  measure.  For  this  change  of  coordinates,  we  have  L0  = 
A  Then  we  consider  the  problem: 

(2.12)  Aw-f  for  y  £  E“  ■  R""1  x  R+ 

Bgw  *  0  for  y  €  R"*1  x  (0). 

Here,  for  simplicity,  we  assume  that  n  >  3  and  introduce  the  notations  y*  = 


'.V' 
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(Yi . y„.1,*yn)  and  wn  the  volume  of  the  unit  ball  in  Rn.  Then,  from  Gilbarg  and 

Trudinger  ([6],  Chapter  6.7),  we  can  solve  (12)  in  terms  of  the  Green’s  f unction: 


w(x)  = 


G(x,y)f(y)dy 


where  G(x,y)  •=  r(x-y)  -  T(x-y*)  +  E(x-y*)  with  T(z)  «  |z|*'n/n(2-n)u>n  the 
fundamental  solution  of  Laplace's  equation  (n  >  2),  and 


_«D 

E(z)  =  -2  c'0(1‘e)  f(z  +  ens)ds,  en-(0 . 0,1). 

*  0 


Since  0  €:  (0,1),  we  have,  for  z£E”  , 


|E(z)|  <  2  f  |-g-  r(z  +  c  s)|ds«2r(z), 
*0 

|w(x)|  <  f  H(x,y)  |  f(y)|  dy  , 

Jc-n 


where  H(x,y)  «  T(x-y)  +  3r(x-y#).  Then,  as  in  Agnon,  Douglis  and  Nirenberg  [1],  we 
can  extend  these  kernels  to  Rn  as  odd  in  xn,  and  apply  the  Calderon-Zygmund 
theorem,  obtaining 

l|w||%  <  Nj  Bf||*j  • 

L  L 

One  can  do  the  same  after  differentiating  w  twice  and  passing  the  derivatives  to 
the  kernel  G,  obtaining  then 


[•] 


(2.13) 


llwir  2  <  Njiini2  j. 

H  L 


/here  the  constants  do  not  depend  upon  0.  Then,  since  LQv  =  Lv  +  (L0-L)v, 


we  have 


||L0v||2  <  2||Lv||2  +  2  ||(L0-L)v||J 

L  L  L 


and  again,  by  (10), 


|L0v||*2  <  2  ||Lv||*  +  6  '||v||2  +  N  '||v||2  2  . 

L  L  H  L 


Thus,  from  (2.13),  considering  a  partition  of  unity  sufficiently  small  so  that 


6'  t  ^Nj)"1,  we  obtain 


||L0v||2  <  2||Lv||2  +6'N2||L0v|!lJ  +  N'||v||*2, 

L  L  L  L 


I|L0v||z2  <  k2(!|Lv||2  +||v||2). 
L  L  L 


So,  again  by  (2.13),  we  sec  that 


i|PiU||22  <  COIPjUll2  2  +  HLpjUll2  2). 
H  L  L 


As  before,  we  can  now  obtain 


llPjUil2  <  K(||g||2  2  +  ||u||J2  +  €||u||*2) 

H  L  L  H 


Pa*  k*  l|*| 

S*!w>»;yV7.»!V 


& 
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3.  An  example. 

In  this  section,  we  discuss  the  situation  in  which  u  in  (1.1),  (1.2)  is  a  scalar, 
fl  =  (0,1)  and  f(u)  is  a  cubic.  It  is  convenient  in  the  computations  to  replace  the 
parameters  (d,0)  by  (dJ,0/(l-0)d).  The  example  to  be  considered  is 

ut  ■  d*uxx  +  f»(u)  *€(0,1) 

(3.1) ,  (l-0)dux  -  0u  at  x  -  0 

(l-0)dux  ■  -0u  at  x  »  1 

where  d  G  (0,°»)  and  0G  [0,1]  and 

(3.2)  f,(u)  -  u(l-u)(u-a),  a  G  [-1,1] 

Since  (3.1)  is  a  gradient  system  every  solution  approaches  an  equilibrium 
solution.  In  the  rescaled  variables  x  -  dy,  these  solutions  satisfy 

(3.3) ,  u^  +  f,(u)  -  0,  y  E  (O.d'1) 

and  boundary  conditions: 


(l-0)uy  -  0u  at  y  -  0 

(3.  4) 

(l*0)uy  -  -0u  at  y  -  d*1. 


Since  the  set  of  equilibrium  solutions  is  bounded,  there  is  a  compact  attractor 
Sd  0  for  every  d  >  0,  0€  [0,1]  (see,  for  example  Henry  [10]  or  Hale  [9]). 

If  Wu(0),  W‘(0)  are  the  stable  and  unstable  sets  for  an  equilibrium  solution 


then  a  recent  result  of  Henry  [11]  shows  that  W"(0)  is  transversal  to  W*(0)  for  all 


equilibrium  solutions  0,0.  This  implies  that  the  flow  defined  by  (3.1)  is 


structurally  stable  if  and  only  if  the  equilibrium  solutions  are  hyperbolic;  that  is, 
if  and  only  if  each  equilibrium  solution  has  the  property  that  its  linear  variational 
equation  has  nonzero  eigenvalues.  This  implies  that  the  curves  in  the  (d,0)-plane 
which  correspond  to  bifurcation  points  of  the  flow  must  be  either  primary 
bifurcations  from  an  equilibrium  or  saddle-node  bifurcations  of  equilibria.  The 
purpose  of  this  section  is  to  discuss  these  curves  for  (3.3),  for  values  of  a  €  [-1,1].  For 
a  =  -1  we  prove  the  following  result  for  the  case  a  ■  -1;  that  is  f(u)  *  u  -  u8. 

Theorem  3.1.  Let  s^  C  (0,«)  x  [0,1]  be  the  structurally  stable  regions  for  (3.1)^ 
which  consists  of  exactly  2j+l  hyperbolic  equilibrium  points.  Then  the  following 
relations  hold: 

1)  Sj  has  only  one  connected  component. 

2)  S0,  Sj  are  unbounded.  Sj  is . t>OV.ndCd.I.pr  j  >  2, 

3)  S0  o  {0  «=  0}  *  <(>,  S0  n  {0  -  1}  *  * 

4)  sJ  n  {0  -  0)  *  4>,  Sj  n  {0  -  l)  *  ft,  Vj  >  l, 
and.  , for  each  imcRgr  k, 

(Ci  uj>k+i  sj>  o  {9  -  0}  -  (Ci  UJ>k  sp  n  {0  -  l)  -  [0,dk], 
where  dk  -  (kit)'1. 

5)  dSj  are  smooth  C^curves  nonincreasing  in  0.  These  curves  arc 
nonintersecting  in  (0,*»)  x  (0,1J. 

Before  proving  this  result,  we  make  the  following  remarks. 


Remark  3.2.  Suppose  dx  is  as  in  property  4)  and  d0  >  dx  is  fixed.  From 
properties  2)  and  5),  if  one  studies  the  attractor  B6q  as  a  function  of  0  for  a  fixed 


d  >  d0,  then  one  must  go  from  a  situation  of  three  equilibrium  points  at  8  ■  0  to 
one  equilibrium  point  at  0  -  1.  Furthermore,  according  to  property  5),  there  is 
only  one  point  0  at  which  a  bifurcation  occurs.  This  d0  provides  a  good  estimate 
of  the  d0  occurring  in  Theorem  1.1. 

Remark  3.3.  Properties  1)  and  4)  imply  that  one  can  find  a  homotopy  from 
any  structurally  stable  system  with  Neumann  conditions  to  a  structurally  stable 
one  with  Dirichlet  conditions.  The  case  with  a  E  (0,i/j)  was  considered  by 
Gardner  [5]  where  he  shows  the  existence  of  such  a  homotopy  for  the  case  with 

three  equilibria.  We  will  see  later  that  for  a  [i/2,lj  no  such  homotopy  exists. 

Proof  of  Theorem  3.1.  Let  L±  be  the  lines  in  the  (u,uy)-plane  defined  by  Lt 

-  {(u,uy)  :  (l-0)uy  -  ±0u).  Let  u  ■  u(y,u„)  be  a  solution  of  (3.3)_1»  where  u0 

corresponds  to  the  maximum  value  of  u.  If  this  maximum  occurs  at  y  ■  T, 
then  utT.Ujj)  «  u0  and  uy(T,u0)  -  0.  We  define  the  "time  map"  T  to  be  T(u0)  «= 
r.  From  this  time  map,  the  existence  of  solutions  of  (3.3)^  can  be  inferred. 
Such  a  solution  exists  if  and  only  if  there  exists  a  u„  E  (0,1)  for  which  T(u0) 

-  (2d)'1.  Introducing  the  polar  coordinates  u  ■  r  cos  s,  uy  »  -r  sin  s  in  (3.3)_r 
one  can  show  that  s  -  s(y,u0)  satisfies  the  differential  equation: 

sy  ■  sin’s  +  (1  -  r’cos*s)cos*s,  y  G  [0,d_1J, 

d«f 

s(0)  -  ♦  -  arctg  0/(l-@)  (0,n/2). 


From  this  equation  we  determine  the  following  expression  for  the  time  map: 

ft 

(3.5)  T^q)  -  [sin’s  +  (1  -  r’cos’sJcosSj^ds 
^  o 
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where  r  «  Ks.Uq).  Then,  as  in  Hale  and  Rocha  [7],  we  can  prove  that  the  time  map 
is  a  monotone  increasing  function  of  u0  (0,1).  In  fact,  differentiating  (3.5),  we 
have: 

(3.6)  T'(Uq)  “  2  f  [sinJs  +  (1  -  racos*s)cosss]'a  r  &  cos4s  ds 

J  o 

and  dr/duQ  >  0  and  t  (0,/r/2)  imply  that  T^Uq)  >  0  for  u0  (0,1).  Thus,  as  in 
Chafee  and  Infante  [3],  the  bifurcation  of  equilibrium  solutions  can  only  occur  at 
the  origin.  This  will  happen  for  the  values  of  d  ■  d(6)  corresponding  to  the  zero 
eigenvalue  for  the  linearized  problem: 

(3.7)  u^  +  u  -  0,  y  (O.d*1), 

and  boundary  conditions  (3.4).  Then,  we  will  have  u  «  A  cos  y  +  B  sin  y,  and  the 
boundary  conditions  will  be  satisfied  if  and  only  if: 

(3.8)  d-  [>«  cotgl^-ji-]] 

This  provides  an  expression  for  the  curves  dSj  referred  to  in  property  5)  of  the 
Theorem  3.1.  Moreover,  one  also  concludes  that  the  dk  in  property  4)  are  given  by 
dk  «  (kn)~*.  This  completes  the  proof  of  the  theorem. 

Figure  3.a  presents  the  sets  Sj,  and  Figure  3.6  the  bifurcation  diagram  for  a 


fixed  value  of  0  (0,1). 


In  the  following,  we  consider  the  time  map  and  its  derivatives  defined  at 
u0  *  0  by  continuity.  Then,  one  can  easily  check  that  T(0)  -  $,  T'(0)  *  0  and 
T"(0)  >  0  for  8  G  (0,1].  For  a  *  -1  the  problem  becomes  much  more  difficult 
because  the  expressions  for  the  derivatives  of  the  time  map  are  more  complicated. 
Nevertheless,  we  can  prove  the  following: 

Theorem  3.4,  For  every  D  sufficiently  large,  let  Sj  c  (0,D]  *  [0,1]  denote., the 
structurally  stable  regions  for  (3.1),  consisting— fi£  exactly  2j+l  h.YP.erh.Q.l?c 
equilibrium  points.  Then  there  exists  a  c  G  (0,1)  such  that  for  all  a  G  (-1,-1+c)  th£ 
following  hold: 

1)  Sj  has  only  one  connected  component  if  j  -  2k,  k  >  0. 

2)  Sj  has  exactly  two  connected  components  if  j  -  2k+I,  k  >  0. 

Moreover,  the  relations  3)  to  3)  of  theorem  3.1  still  hold. 

To  prove  this,  we  introduce  in  the  time  map  the  dependence  on  a,  T  «  T(u„,a): 

f  ^ 

(3.9)  T(u„,a)  -  (sin*s  +  [-a  +  (l+a)r  cos  -  r,cos*s]cos*s)‘1ds 

From  the  remark  before  the  statement  of  the  theorem,  we  know  that  0T/9uo(O,-l)  - 
0  and  9*T/9u2(0,-1)  >  0  for  8  G  (0,1  j.  In  the  same  way,  one  can  verify  that 
O^T/daOu^O.-l)  <  0  for  8G  (0,1].  Then,  for  any  8  >  0,  we  can  find  an  c  >  0  such 
that,  for  8  [8,1],  we  have  ^T/Ou^O.-l)  >  «  and  0*T/9a8uo(O,-l)  <  Hence,  for 

6  G  [8,1],  the  changes  introduced  in  the  time  map  as  a  >  *1  are  very  simple,  and 
we  can  find  a  c  G  (0,1]  such  that  for  all  a  G  (-1,-1 +c)  the  time  map  has  a  unique 
extremum  at  u0G  (0,1),  which  is  a  minimum.  This  gives  us  the  shape  of  the  first 
bifurcation  curve  in  Figure  4.a,  showing  what  is  usually  called  a  transcritical 
bifurcation  at  the  origin.  A  simple  analysis  of  the  phase  plane  shows  that  only 
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the  odd  bifurcations  at  the  origin  will  be  transcritical,  the  even  ones  being 
supercritical.  This  observation  takes  care  of  the  curves  9S^  in  the  region  [0,D]  * 
(6,1].  For  the  region  (0,D)  x  (8,1],  we  start  by  observing  that  if  8  -  6(D)  is  small 
enough  this  region  always  contains  at  least  three  hyperbolic  equilibria,  thus,  the 
first  bifurcation  of  the  origin  is  excluded.  Then,  one  needs  only  to  consider  the 
solutions  arising  from  the  second,  third,  etc.,  bifurcations.  If  we  define 
d«f  f  *+U-1)n 

Uj(u0,a)  -  (sin*s  +  [-a  +  (l+a)r  cos  s  -  r*cos,a]cos*s)’1ds, 

* 


1  2 
1  »*»••• 


then  Uj  represents  the  value  of  y  at  which  the  solution  of  (3.3)a  satisfying  the 
initial  condition  in  (3.4)  and  having  at  the  first  maximum  the  value  u0,  satisfies 
the  final  condition  in  (3.4)  after  passing  through  j  extrema.  Note  the  relation 
with  the  time  map:  U^Ug.a)  «  IT^a).  One  can  clearly  use  Uj  to  determine  the 
existence  of  solutions  of  (3.3)a  in  the  same  way  as  the  time  map  was  used.  As 
before,  we  can  verify  now  that,  for  all  8  E  [0,1],  9*Uj/9uJ(0,-l)  >  0  for  j  >  2, 
3Uj/9uo(0,a)  -  0  for  j  -  2k  and  all  a  >  -1,  and  3*Uj/3la3uo(0,-l)  <  0  for  j  -  2k+l,  k 
-  1,2,...  .  Hence,  the  changes  introduced  in  Uj,  j  >  2,  as  a  >  -1  are  very  simple  and 
again  we  can  find  a  c  E  (0,1]  such  that,  for  all  a  E  (-1,-1+c),  Ujt  j  >  2,  has  a 
unique  extremum  which  is  a  minimum.  This  minimum  occurs  at  the  origin  if  j  * 
2k,  and  at  Uj  E  (0,1)  if  j  -  2k+l,  for  k  -  1,2,...  .  This  justifies  the  bifurcation 
diagram  presented  in  Figure  4.a,  and  concludes  the  proof  of  the  theorem.  In 
Figure  4.b,  we  present  the  sets  Sj  as  obtained  from  the  theorem.  It  turns  out  that, 
if  we  consider  the  linearized  problem  -  au  »  0  and  compare  with  (3.7),  we 
obtain  an  expression  for  the  curves  3Sj  corresponding  to  the  bifurcations  at  the 
origin  if  we  multiply  (3.8)  by  the  factor  (-a)1/*  for  a£[-l,0). 


WV,1 


-28- 


Figure  4.a  Figure  4.b 

If  one  considers  the  results  obtained  by  Smoller  and  Wasserman  [IS]  for  the  cases 
of  Dirichlet  and  Neumann  boundary  conditions,  the  results  of  this  theorem  are  not 
surprising.  Moreover,  numerical  tests  indicate  that  these  results  seem  to  hold  for 
all  a  in  [-1,0);  thus,  the  maximum  value  of  c  in  theorem  3.4  being  possibly  1. 

For  a  -  0,  the  problem  is  degenerate  and  does  not  have  any  structurally 
stable  regions.  It  becomes  then  very  interesting  to  make  the  same  study  for  a£ 
(0,1).  This  problem  is  as  difficult  as  the  previous  one  for  a  £  [-1.0)  for  the  same 
reasons.  Therefore,  we  concentrate  on  qualitative  information.  Considering  the 
phase  diagram  corresponding  to  the  equation  (3.3),,  one  notices  that  there  is  a 
qualitative  change  as  a  crosses  the  value  1/2.  In  fact,  for  a  £(0,1/2),  this  diagram 
contains  a  homoclinic  orbit  to  the  point  (u,uy)  -  (0,0);  at  a  -  1/2,  it  contains  two 
heteroclinic  orbits  to  the  points  (0,0)  and  (1,0),  and,  for  a  £  (1/2,1),  the  diagram 


has  an  orbit  homoclinic  to  the  point  (1,0).  This  qualitatiave  change  reflects  on  the 
structurally  stable  regions  for  (3.1),.  Let  us  define  a  function  7  «  7(a)  in  the 
following  way.  For  a  G  (0,1/2],  7(a)  is  the  value  of  0  in  L±  corresponding  to  the 
angle  of  the  tangents  to  the  separatrices  at  the  origin  in  the  phase  diagram.  A 
simple  computation  yields  7(a)  -  (1  +  a"1/*)*1-  For  a  G  (1/2,1)  we  define  7(a)  as 
the  value  of  8  in  L±  corresponding  to  the  tangents  to  the  homoclinic  orbit  passing 
through  u  -  1.  This  function  7  is  a  continuous  function  satisfying  0  <  7(a)  < 
(l+ziT)*1  for  a  G  (0,1).  Then,  if  again  we  let  Sj  denote  the  structurally  stable 
regions  for  (3.1),  corresponding  to  2j+l  equilibria,  we  have: 

Theorem  3.5:  For  a  G  (0.1)  the  following  holds: 

1)  If  a  <  1/2.  there  exist  positive  constants  6  and  D  such  that 
(0,8)  x  [7(a),  1]  C  Sv 

ID,-)  x  [7(a),  1]  C  Sq. 

2)  1£  a  >  1/2,  wfi  Myg 
(0,-)  x  [7(a),  1]  C  S0. 

To  prove  this,  we  consider  the  effect  that  the  above  observation  about  the 
phase  diagram  for  (3.1),  has  upon  the  time  map  as  defined  by  (3.9).  If,  for  a 
(0,1/2),  we  denote  by  a  the  smallest  positive  root  of  Jj}  f,(s)ds  -  0,  then  the  point 
(ce,0)  of  the  phase  diagram  corresponds  to  the  intersection  of  the  homoclinic 
orbit  through  the  origin  with  the  positive  u-axis.  Then,  for  6  )  7(a),  the  time  map 
T(-,a)  :  (a,l)  -»  (0,-)  is  continuously  differentiable  and  unbounded  as  u0  -*  o  or  1. 
Moreover,  one  can  easily  check  that  3T/9u0(-,a)  :  («,!)  •*  R  is  positive  as  u0 
approaches  1  and  negative  as  u0  approaches  a  This  implies  the  existence  of  a 
minimum  value  p0  >  0  for  T(-,a),  and  also  a  maximum  value  p,  a  p0  for  its 
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Remark  3.6.  It  is  clear  that  part  2)  of  the  previous  theorem  presents  the 
existence  of  a  homotopy  from  a  structurally  stable  system  with  Neumann 
boundary  conditions,  which  must  have  at  least  three  solutions,  to  a  structurally 
stable  system  with  Dirichlet  boundary  conidtions,  which  must  have  only  one 
solution.  The  above  example  also  makes  clear  what  one  should  do  to  create  such 
qualitative  phenomena  as  the  alternative  in  Theorem  3.5,  in  systems  corresponding 
to  more  general  functions  f.  Finally,  in  this  example  the  curves  cK.  in  theorems 
3.1  and  3.4  were  constituted  only  by  codimension  one  bifurcations.  It  is  possible 
to  create  examples  in  which  these  lines  intersect,  presenting  higher  codimension 
bifurcations.  For  instance,  take  an  example  for  which  the  time  map  at  some  90 
has  two  minima  which  are  equal.  Then,  for  this  example  there  would  be  a  d0  >  0 
such  that  (do,0o)  would  correspond  to  a  codimension  two  bifurcation. 


